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ABSTRACT 

We assessed the effects of deforestation on soil carbon (C) and nutrient stocks in the premontane landscape near Las Cruces Biological 
Station in southern Costa Rica, where forests were cleared for pasture in the mid-1960s. We excavated six soil pits to a depth of 1 m in 
both pasture and primary forest, and found that C stocks were —20 kg C/m“ in both settings. Nevertheless, soil <5 13 C suggests —50 per- 
cent of the forest-derived soil C above 40 cm depth has turned over since deforestation. Soil nitrogen (N) and phosphorus (P) stocks 
derived from the soil pits were not significantly different between land uses (P = 0.43 and 0.61, respectively). At a larger spatial scale, 
however, the ubiquity of ruts produced by catde-induced erosion indicates that there are substantial soil effects of grazing in this steep 
landscape. Ruts averaged 13 cm deep and covered —45 percent of the landscape, and thus are evidence of the removal of 0.7 Mg C/ 
ha/yr, and 70, 9 and 40 kg/ha/yr of N, P and potassium (K), respectively. Subsoils in this region are —10 times less C- and N-rich, and 
—2 times less P- and K-rich than the topsoil. Thus, rapid topsoil loss may lead to a decline in pasture productivity in the coming dec- 
ades. These data also suggest that the soil C footprint of deforestation in this landscape may be determined by the fate of soil C as it is 
transported downstream, rather than C turnover in situ. 

Abstract in Spanish is available in the online version of this article 
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Tropical forests are being cleared at unprecedented rates, 
driven largely by the demand for food from an increasing popu- 
lation, changing diets, and economic development (Smith et al. 
2007, DeFries et al. 2010). Tropical deforestation has accounted 
for —60 percent of global emissions due to land-use change since 
1850 (Houghton 2003), and today constitutes —15 percent of cur- 
rent anthropogenic emissions (Houghton 2007). In Central and 
South America, deforestation is largely driven by pasture and 
agricultural expansion, and in many regions these practices oper- 
ate with minimal (if any) fertilizer addition. With low inputs, land 
degradation post-deforestation can be rapid (McGrath et al. 
2001), and as nutrient losses come to constrain production, there 
is an incentive to clear more land (Serrao et al. 1996). Within this 
context, however, there is considerable variation in outcome given 
the large differences in soil fertility across the tropics (Sanchez 
1976). 

The replacement of productive forests by agriculture or pas- 
ture results in substantial carbon (C) emissions from above- 
ground biomass. In contrast, recent meta-analyses find relatively 
small effects on soil C ( eg ., Powers et al. 201 1), with large vari- 
ance around the mean effect of a particular transition {e.g., forest 
to pasture). High spatial heterogeneity in both rainfall and soil 
type, and clustering of studies in particular regions, preclude a 
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robust extrapolation of these results to the tropics as a whole 
(Powers et al. 2011). Furthermore, many studies report data 
from the shallow soil (often <30 cm depth e.g., Schedlbauer & 
Kavanagh 2008), which may not capture the integrated effects 
throughout the soil column (Veldkamp 1994, Guo & Gifford 
2002). Perhaps more concerning, most studies focus on changes 
in carbon concentration, which is only relevant to global change 
if bulk density is measured accurately (which is uncommon); 
otherwise compaction in pastures may lead to a substantial over- 
estimate of soil carbon content (Veldkamp et al. 2003). Finally, it 
is important to note, however, that even pools of soil C with no 
net change may be turning over, as demonstrated by heavier 
stable isotope signatures (<5 13 C) in pasture soils where grasses 
utilizing C4 photosynthesis have replaced C3 forest vegetation 
(Veldkamp 1994). Thus, our understanding of changes in soil C 
stocks associated with forest conversion to pasture bears further 
consideration, and is one of the foci of this article. 

Land-use change in the tropics can also have significant 
effects on soil fertility and aboveground productivity (Davidson 
et al. 2004), and many tropical grazing operations run measureable 
nutrient deficits (Smaling et al. 1997, MacDonald et al. 201 1). 
Nutrient budgets, however, are difficult to scale up to the 
landscape (Fearnside 1980, Smaling et al. 1997, Hartemink, 2006), 
in part because of poorly constrained estimates of soil erosion 
(Davidson et al. 2004, Lai 2009). Even where erosion has been 
considered, the majority of studies have focused on relatively 
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stable tectonic settings with low inherent soil fertility and very low 
erosion rates, such as sub-Saharan Africa and the Amazon basin 
(McGrath et al. 2001, Sanchez 2002). Studies of nutrient stock 
changes associated with deforestation for pasture are rare in the 
montane tropics, where erosion rates are likely to be higher. 

In this context we measured changes in total soil carbon 
(C), nitrogen (N), phosphorus (P), and potassium (K) associated 
with deforestation for pasture in the premontane forests near Las 
Cruces Biological Station in Southern Costa Rica. We worked at 
two spatial scales. First, we excavated multiple soil pits and we 
measured soil C, N, P, and K content in a primary forest and 
pasture cleared —50 yr ago. We used <5 13 C of the soil and forest 
and pasture vegetation to assess how much of the soil C stock 
has turned over since deforestation. At a larger spatial scale, we 
quantified topsoil erosion and C, N, P, and K exported from 
a network of extensive cattle ruts that cut across these steep 
pastures. 

METHODS 

Site description. — We worked in and around Las Cruces Biolog- 
ical Station in Coto Brus County, southern Costa Rica in 2009— 
2010 (Fig. SI). The station has a mean annual temperature of 
21°C, and mean annual rainfall is —4000 mm/yr, with a pro- 
nounced dry season (<100 mm/mo) from December to April. 
The study area is in the foothills of the Zapote range at an alti- 
tude of 1200—1300 m, and is dominated by steep topography 
(>20° in many places). Soils are mostly Lithic and Andie Dys- 
trudepts along with Typic Hapludands (Jin et al. 2000) derived 
from basaltic andesite and associated volcanics. The region 
hosted premontane tropical forest until the mid-20th century, 
when rapid deforestation began, and it is now dominated by pas- 
ture and small agricultural holdings. Primary forest remains in 
patches, of which Las Cruces is one. 

We selected sites adjacent to Las Cruces Biological Station 
that were converted direedy from forest to pasture in the 1960s 
and remained actively grazed pasture until the station purchased 
the 26 ha property in 2008 (Z. Zahawi, pers. comm.). Vegetation 
was not burned after clearing and we found no evidence of char- 
coal while sampling soils. These pastures were seeded soon after 
deforestation, and such ‘low-investment systems’ are typically not 
fertilized (Jin et al. 2000). The pastures are now dominated by 
Axonopus micay and Urochloa mollis, both C4 grasses of Neotropical 
origin. 

Analyses: Soils. — In August 2009, we excavated 12 soil pits, six 
in pasture and six in forest, arrayed across —600 m of both land 
use types. Pasture pits were located on a single 26 ha parcel of 
land adjacent to the forest. We sampled on three separate ridges 
and slopes, and all pasture pits were within 200 m of the forest 
edge. Likewise, all forest pits were located within one continuous 
parcel of forest (Las Cruces Biological Station), included ridge 
(N = 3) and slope (N = 3) positions, and were between 250 and 
500 m from the pasture edge. We dug pits to a depth of 1 m. 
Once each pit was open, we excavated a 10 x 10 cm shaft down 


the sampling face to assess bulk density, and sampled the soils 
from this shaft for additional chemical analyses. The soils showed 
strong horizonation so we chose to sample by horizon rather 
than depth. Where horizons were >20 cm thick we took several 
within-horizon subsamples at different depths and analyzed them 
separately, we report means by horizon below. We air-dried sam- 
ples at Las Cruces and shipped homogenized subsamples back to 
Brown University, where they were sieved (<2 mm) and hand- 
picked for roots. We split each sample into three subsamples. 
One was analyzed for total elements by lithium borate fusion and 
ICP-MS (ALS Chemex, Reno, Nevada, U.S.A.), a second for C, 
N, and <5 13 C on a Costech Analytical Model 4010 (Costech Ana- 
lytical Technologies, Valencia, California, U.S.A.) flash-combustion 
element analyzer with a Thermo Conflo III interface into a Del- 
ta VPlus isotope ratio mass spectrometer (Thermo Scientific, West 
Palm Beach, Florida, U.S.A.) in the Whiteside lab at Brown Uni- 
versity, and a third for bulk density by drying at 105°C for 48 hr. 
All results are presented on an oven dry weight basis. 

We also collected leaves and litter from both pasture and 
forest, which were dried at 65°C and analyzed for C, N, and 
<5 13 C as described above. We calculated the contribution of C3 
(forest) and C4 (pasture grasses) plants to soil carbon as: 


^pastureC ' 
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where F pas tureC is the fraction of soil C from pasture inputs (Pow- 
ers & Veldkamp 2005). There are several complications to this 
approach discussed in detail by Powers and Veldkamp (2005), 
including the need to measure both leaf and root litter, and the 
presence of both C3 and C4 plants in pasture. We assumed no 
C3 inputs to the pasture after its creation, and did not measure 
the <5 13 C of roots in either forest or pasture, so our calculation of 
Fpasturec may be less accurate than one that takes these factors 
into consideration. We did correct for changes, however, in 5 ’C 
of forest soils with depth in our calculation of F pasture c. We esti- 
mated S u Ci OICSt so fl for depth ( D ) in pasture from a third order 
polynomial derived from the observed relationship between 
d 13 C (oast soil and D at our six forest soil pits (R“ = 0.74): 

<5 13 C d = -25.9 + 0.0465D - 0.000812D 2 + (3.59exp“ 6 )D 3 (2) 

We then used as the forest soil end member for each 

depth at which we calculated F pastureC . 

Analyses: Soil loss from cattle ruts. — Our initial observations 
of the landscape indicated a high density of ruts in pasture on 
the slopes >20° that were absent from forest slopes (Fig. S2). 
Such features have been noted for almost a century (Odum 1922 
quoted in Rahm 1962), and are now generally considered animal 
paths dubbed ‘terracettes’ or ‘cowtours’ in the geomorphic litera- 
ture (Trimble & Mendel 1995). On the pasture hillslopes near 
Las Cruces, ruts run approximately perpendicular to slopes, and 
the uphill edge of one rut is often the downhill edge of the rut 
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above it, forming a hiUslope that resembles a flight of stairs. We 
measured the dimensions of individual ruts (N = 87) down seven 
slopes, as well as the fraction of the hillslopes covered by ruts. 
Most ruts could be approximated as triangular cross-section cuts 
into the hillslope (Fig. S2), although occasionally, ruts were carved 
deeper, resulting in trapezoidal cross-sections. We measured the 
dimensions of the ruts, and calculated their mean depth relative 
to the original hillslope. This allowed us to estimate total soil loss 
associated with the ruts. We converted depth of soil loss to mass 
loss (kg/m - ) by multiplying depth by the mean bulk density of 
forest soil for a given depth. 

To test the hypothesis that the ruts were generated by ero- 
sion rather than compaction, we took soil cores from the upslope 
and downslope edge of several ruts (0—5, 5—10, 10—15 cm depth 
below the rut surface). As soil percent carbon (%C) decreases 
with depth in almost all soils Qobbagy & Jackson 2001), we 
tested our hypothesis by comparing soil %C between the soils 
collected from the uphill and downhill edges of the ruts. If 
erosion was the driver of rut formation, we expected that soil at 
the uphill edge of ruts would exhibit characteristics of deeper soil 
( e.g. ., lower C concentration) than cores from downhill edge 
because more soil had been removed above the exposed soil near 
the uphill edge of the rut. In contrast, if the ruts were formed by 
compaction, we expected no difference in soil C concentration 
between the uphill and downhill edges of the rut. 

As a second test of the hypothesis that these ruts were 
caused by erosion, rather than compaction, we used our calcula- 
tion of rut depth (based on the dimensions of the rut) to gener- 
ate an ‘original’ depth of each rut soil sample. We define original 
depth as the depth of the soil before the formation of the rut 
{e.g., if the rut is 10 cm deep at the point of sampling, and we 
sampled 0-5 cm below the surface, we assigned an ‘original’ 
depth of 12.5 cm). We then used the observed relationship 
between soil depth and soil %C from our soil pits (Fig. 1A) to 
generate an ‘expected’ %C for the rut soils given their calculated 
original depth in the soil profile. This expected soil %C should 
only correlate with measured C if the ruts were formed by ero- 
sion. Thus, we took a significant correlation between expected 
and measured soil %C as evidence in support of this hypothesis. 

Statistics. — Carbon, N, and <5 13 C between forest and pasture 
were compared using restricted maximum likelihood estimation 
with land use type, depth and landscape positions as fixed effects. 
Comparisons between the soil %C in the front and back of ruts 
was done using one-way AN OVA after testing data for normality 
(Shapario— Wilk W test; P > 0.05 for all depths on both the uphill 
and downhill edges of the ruts). To calculate the expected %C 
for a soil in a rut of a given depth (D), we fitted a first order 
decay function to the relationship between soil %C and D 
(R 2 = 0.89) such that: 

%C = 8.49x«v>(-0.042xD) + 0.8 (3) 

We then calculated an expected %C based on the depth of 
soil sampled in a particular rut plus the soil we calculated had 
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FIGURE 1. (A) Soil C concentrations; and (B) rl 1 'Cl. Filled squares are forest 
pits, open triangles are pasture pits. 


been removed based on the geometry of the rut. Expected and 
measured %C were compared by Pearson’s correlation. All statis- 
tical analyses were performed in JMP (v. 8.0.2, SAS Institute, Inc. 
2009; Cary, North Carolina, U.S.A.). 

RESULTS 

Soil carbon. — Carbon concentrations decreased from —10 per- 
cent near the surface to —1 percent at a depth of 1 m, but did 
not differ between forest and pasture (Fig. 1A). Neither did bulk 
density (One-way ANOVA, df =10, P = 0.64; Table SI), and 
thus total soil C stocks to a depth of 1 m were similar between 
forest and pasture (P = 0.79; mean —20 ± 0.8 kg C/m - for both 
land uses). Soil within the top 10 cm contained about 6.0 kg C/ 
m - , and soil between 90—100 cm depth contained 0.60 kg C/m - . 
Although the trends in soil C concentrations are consistent across 
all pits, variation in bulk density makes the total C content quite 
variable across the landscape (Table SI). 

Forest c) 13 C values increased from about — 27% 0 near the 
surface to between — 25% 0 and — 21% 0 at 1 m depth (Fig. IB). In 
pasture, 5 C of bulk soil decreased from about — 17% 0 at the 
surface to values indistinguishable from forest at 1 m depth. 
<5 13 C in the A horizon differed significantly between forest and 
pasture (P < 0.0001), but (5 13 C of deeper soils (B horizons) in 
both settings were not significantly different (P = 0.19), converg- 
ing to a value of about — 23% 0 . Litter from the C3 trees in the 
forest has a <5 L ’C — 30% o , whereas the C4 grass in pastures has 
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values ~12% 0 . From equations (1) and (2), the uppermost pasture 
soils have derived —50-60 percent of their carbon from C4 
sources since deforestation, whereas at 1 m depth no C4 carbon 
is detectable. 

Soil nutrient content. — Soil N decreased with depth, from 
between 0.6 and 1 percent of total soil mass at the surface to about 
0.1 percent at a depth of 1 m (Table SI). Soil N stocks decreased 
with depth from about 0.60 kg N/m 2 within the top 10 cm to 
0.042 kg N/m" between 90 and 100 cm depth (Table SI). Total 
soil N stocks to 1 m did not differ between forest and pasture 
(1.9 ± 0.1 and 1.7 ± 0.2 kg N/m - , respectively, P = 0.43). Soil P 
contents decreased from 800 to 1700 ppm at the soil surface to 
about 200—1000 ppm at 1 m, although depth profiles were more 
variable than either C or N (Table SI). Mass of P varied between 
0.076 kg P/m 2 in the top 10 cm and 0.033 kg P/m 2 at 90— 
100 cm depth. Total P stocks to 1 m were not different between 
forest and pasture (0.57 ± 0.10 and 0.55 ± 0.08 kg P/m 2 , respec- 
tively, P = 0.89). Soil K decreased from 0.3 to 0.7 percent at the 
surface to about 0.15 percent at 1 m (Table SI). Total soil K stocks 
to 1 m did not differ between forest and pasture (2.7 ± 0.5 and 
2.1 ± 0.1 kg K/m 2 , respectively, P = 0.35). The full soil dataset is 
available in the supplemental online material. 

Soil loss from cattle ruts. — Although our analyses at the scale 
of a soil shaft (10 cm x 10 cm x 1 m) showed no differences 
between forest and pasture C and nutrient stocks, the presence 
of extensive ruts in the pasture suggested that the scale of a soil 
pit was too small to capture the relevant losses of C and nutri- 
ents. Ruts measured 13 ± 0.5 cm deep, and over the seven mea- 
sured hillslope transects extended over 45 percent of the 
landscape. Our observation of red B-horizon soils in the bottom 
of many ruts led us to postulate erosion as the rut-forming 
agent, which generated the expectation that soils in the uphill 
portion of the rut should be more similar to soils from lower in 
an undisturbed profile than soils from the downhill edge of the 
ruts. Consistent with this hypothesis, soil C in the 15 cm of soil 
below the rut surface was significantly lower in soils from the 
uphill portion of the rut (one-way ANOVA; df = 46, 
P = 0.006). Although this was not true of the upper 5 cm of 
rut soil (Tukey’s HSD, df = 14, P > 0.05) it was true for both 5 
—10 and 10—15 cm depth (Tukey’s HSD, df = 14, P = 0.02 in 
both cases). 

Using the measured geometry of the ruts, we generated an 
‘original’ soil depth for each rut core, and used equation (3) to 
translate this value into an expected soil %C based on the empiri- 
cal relationship between forest soil C and depth (Fig. 1A). 
Expected soil C was significandy correlated with measured soil % 
C in the ruts ( N = 48, r = 0.46, P = 0.0005; Fig. 2), further sup- 
porting the hypothesis that erosion, rather than compaction, had 
formed these prominent features in the landscape. When analyz- 
ing each depth interval separately, the 0-5 cm rut soils showed 
no significant correlation between predicted and observed soil C 
(df = 14, r = 0.09, P = 0.7), whereas both the 5-10 and 10— 
15 cm intervals showed a strong correlation between expected 
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FIGURE 2. Measured soil %C from the front and back of ruts plotted 
against the predicted soil %C based on the relationship observed between soil 
%C and depth in our soil pits. The overall relationship is significant (r = 0.48, 
P= 0.0005), but is driven by soil from 5 to 10 (open triangles, r = 0.61, 
P = 0.01) and 10 to 15 cm (gray circles, r = 0.58, P = 0.02) beneath the rut 
surface. The relationship in the 0—5 cm rut soil is not significant (black 
squares, r— 0.1, P = 0.7). 


and measured soil carbon (r = 0.62, 0.58, P= 0.01, 0.02, respec- 
tively, df = 14 for both). 

DISCUSSION 

These data suggest that the montane soils of Las Cruces hold 
substantially more carbon in the top 30 cm (—13 kg C/m 2 ) than 
several other lowland tropical sites that have been intensively 
studied. For example, the lowland (<100 m asl) soils near La 
Selva Biological Station in northern Costa Rica store 9.7 and 
5.4 kg C/m 2 (Veldkamp 1994), and a broader exploration of 
forests and pastures in the same region averaged —8 kg C/m 2 — 
30 cm depth (Powers & Veldkamp 2005). Our sites are at 
higher elevation than those (although on similar parent mate- 
rial), so it is not surprising that soil C content is higher as well. 
The role of parent material and mineralogy may help explain 
why the soil C content around Las Cruces exceeds that in a 
similar elevation lower montane system in Ecuador, where soils 
contain between 8 and 10 kg C/m 2 to a depth of 30 cm 
(Rhoades et al. 2000). In other, more highly weathered lowlands 
such as the Amazon, C stocks are also considerably lower. A 
study by Veldkamp et al. (2003) found that soil C contents at La 
Selva were twice as high as those reported from the Amazon 
basin, and stocks in Rondonia range from 2.4 to 6.2 kg C/m“ 
in the top 30 cm of forest soils (Neill et al. 1997). Although the 
soils at Las Cruces hold large stocks of C, this may not be 
uncommon for soils at similar elevation and derived from simi- 
lar parent material. 



Soil Carbon and Nutrient Changes 665 


At the scale of a sampling shaft the soil C and nutrient 
stocks at Las Cruces appear unaltered by forest conversion to 
pasture. This result is consistent with a meta-analysis of carbon 
changes associated with deforestation in tropical forests with 
over 3000 mm/yr of precipitation; indeed, many of the data in 
this category come from other studies in premontane Costa Rica 
(Guo & Gifford 2002). Among other studies conducted in 
Costa Rica, Veldkamp (1994) found decreases in total soil C fol- 
lowing conversion of forest to pasture in Atlantic Costa Rica, 
particularly in the top 10 cm, but ‘passive’ C stocks at greater 
depths were more stable. In contrast, several authors have 
observed no significant changes in soil C stocks to 30 cm depth 
in pasture after deforestation across a range of pasture ages in 
northeastern Costa Rica (Powers & Veldkamp 2005, Schedlbauer 
& Kavanagh 2008). Powers and Veldkamp (2005) studied sites 
with similar annual precipitation (—4000 mm) to Las Cruces, but 
which were located at lower elevations (<750 m) and thus had 
warmer temperatures (about 21— 25°C). They noted that despite 
no significant mean change between land uses, individual site 
changes in soil C stocks were variable in both sign and direc- 
tion, as has been found in the Amazon (Neill et al. 1997). The 
lack of change in C stocks at Las Cruces does not preclude C 
turnover, however. <5 13 C suggests that —50 percent of the soil C 
in the upper 40 cm of pasture has turned over since deforesta- 
tion. Furthermore, if the loss of soil nutrients from ruts 
adversely affects pasture productivity, soil C loss may ensue as 
carbon inputs to soils diminish. 

Although the sampling shafts show no effect of deforesta- 
tion on soil C and nutrient stocks, the network of ruts across 
the pasture landscape suggests a mechanism for potentially large 
soil loss at the scale of a hillslope. The resolution of this appar- 
ent contradiction can be found in the scale of our soil pits rela- 
tive to the landscape. When digging soil pits on steep slopes, 
the sampling face is often (and in our case was always) located 
on the uphill section of the pit. Although this minimizes dig- 
ging time, it increases the likelihood that the sampling face is 
cut into the downhill side of a rut (e.g., along the dashed blue 
line in Fig. S2). This profile location will have the least soil loss 
of anywhere in this heavily rutted landscape, thus it gives the 
most complete soil profile available in the pastures. However it 
also leads to an underestimate of the dramatic changes associ- 
ated with grazing in this landscape. In our system, ruts average 
13 cm deep and cover 45 percent of the hillslopes we mea- 
sured, which translates into an erosion rate of —1.3 mm/yr just 
from ruts. For comparison, long-term erosion rates determined 
by cosmogenic radionuclide analyses in forested watersheds at 
nearby Las Alturas station are —0.2 mm/yr (S. Porder, unpubl. 
data). Our data are consistent with erosion as the dominant 
process by which these ruts are created, and the geomorphic lit- 
erature on ‘terracettes’ and ‘cowtours’. We performed two inde- 
pendent tests of the hypothesis that the ruts are created by 
erosion and compaction: (1) comparison of soils from the uphill 
and downhill portion of individual ruts; and (2) a calculation of 
expected rut soil C based on the geometry of the rut, the 
assumption that ruts had been formed by erosion, and compari- 


son to measured soil C. Both tests support the erosion hypoth- 
esis, although the 0—5 cm expected %C data do not, suggesting 
compaction in these upper soils. As the upper 5 cm of soil can 
also be affected by topsoil washing downslope (which would 
increase and homogenize soil C in the very upper soil of the 
rut), we believe the 5—10 and 10—15 cm data are more indica- 
tive of the true differences between the uphill and downhill 
edges of the ruts. 

Our data suggest that —6 cm of soil have, on average, been 
lost from these pastures via erosion since deforestation ca 1960. 
Based on soil nutrient stocks in the upper 6 cm, we estimate that 
since forest conversion 32 Mg C/ha, or —18 percent of the total 
soil C stock to a depth of 1 m has been removed from the site 
and washed downstream. Similarly, —3500 kg N/ha (18% of 1 m 
stock), 450 kg P/ha (9% of 1 m stock), and 2000 kg K/ha 
(7.6% of 1 m stock) have been lost. This translates into a rate of 
C and nutrient loss of 0.64 Mg C/ha/yr, 70 kg N/ha/yr, 
9 kg P/ha/yr, and 40 kg K/ha/yr. 

We suggest that these losses, rather than exports in beef or 
from decomposition, will drive the larger impact of deforestation 
on soil carbon and nutrient stocks in this wet, premontane land- 
scape. Globally, the fate of eroded carbon (burial and storage vs. 
instream processing and C release to the atmosphere) is hotly 
debated (Lai 2003, Lai et al. 2004, Van Oost et al. 2007). Because 
we do not yet know the fate of this eroded soil C at Las Cruces, 
we cannot quantify whether rut-driven erosion in Las Cruces is 
ultimately a source or sink of soil C. 

The importance of erosion in this landscape extends 
beyond the carbon cycle. The rapid loss of fertile topsoil may 
come to constrain pasture productivity around Las Cruces in 
the coming decades, as it has elsewhere (Hartemink 2006). 
Depletion of soil nutrients has significantly impacted global 
crop yields (Tan et al. 2005, Lai 2009), and we expect the 
losses observed at Las Cruces to lower pasture productivity 
given the much lower N and P found in the subsoil. Although 
N and P enrichment of topsoil relative to subsoil is a global 
phenomenon even in unfertilized settings (Jobbagy & Jackson 
2001), N enrichment can occur naturally over the course of 
several decades (Davidson et al. 2007). In contrast, P enrich- 
ment by plant cycling can take tens of thousands of years 
(Porder & Chadwick 2009). The P concentration in these pas- 
tures is 4—5 times what is found in similar settings in the 
Amazon (where pasture productivity degrades over a decade or 
two), so it is no surprise that these pastures have remained 
fertile for over 50 yr. But the pace of cattle-induced soil loss, 
and the lower N and P concentrations in the exposed subsoil, 
suggests that productivity is likely to diminish under this man- 
agement regimen in the coming decades. Although we caution 
against blanket extrapolation of our results given the limited 
spatial scale of our sampling, cattle ruts are ubiquitous in 
almost every steep pasture we have observed in southern 
Costa Rica (and indeed across much of the tropics). Thus, 
we suspect that this phenomenon may be an important one 
to consider in evaluating the soil carbon and nutrient effects 
of deforestation. 
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CONCLUSIONS 

In some ways, these data fit many well established patterns from 
the literature: soil carbon does not change with forest conversion 
to pasture, montane tropical soils are C rich relative to their low- 
land counterparts, and soil nutrient stocks are high on volcanic 
soils. The ubiquity of cattle ruts, however, both in Las Cruces 
and throughout the tropics, suggests that perhaps a focus on car- 
bon and nutrients at the scale of a soil pit or soil core is not 
appropriate in this landscape. At the landscape scale, grazing ani- 
mals are causing dramatic shifts in both C and nutrients, and the 
fate of that C, and of pasture productivity in the longer term, will 
almost certainly be affected by this process. Although this study 
is only a start in quantifying the extent and impact of this larger- 
scale landscape feature, we suspect it may well be important for 
our understanding of the effects of land use on tropical soil 
properties and management. 
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